Exp-function method for solving the Burgers-Fisher equation with variable coefficients 
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In this paper, the exp-function method with the aid of symbolic computational system is used to obtain gen- 
eralized travelling wave solutions of a Burgers-Fisher equation with variable coefficients. It is shown that the 
exp-function method, with the help of symbolic computation, provides a straightforward and powerful mathe- 
matical tool to solve the nonlinear evolution equation with variable coefficients in mathematical physics. 
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I. INTRODUCTION 



The investigation of exact solutions of nonlinear evolution equations(NLEEs) plays an important role in the study of nonlinear 
physics phenomena. The importance of obtaining the exact solutions of these nonlinear equations, if available, will facilitate 
T— I the verification of numerical solvers and aids in the stability analysis of solutions. In the past several decades, many effective 
7—i methods for obtaining exact solutions of NLEEs have been presented, such as the tanh-function method [B-tlll , extended tanh 
method |4, 5], F-expansion method |6, 7], sine-cosine method |8-10], Jacobian elliptic function method j 1 ll - [T4ll . homotopy 
|~~J ^ perturbation method 1 15- 18], variational iteration method 1 19, 20] and Adomian method ]21-23] and so on. 

Recently, He and Wu 124] proposed a straightforward and concise method, called Exp-function method, to obtain generalized 
solitary solutions and periodic solutions. Applications of the Exp-function method can be found in ]25- 27] for solving nonlinear 
Ph evolution equations arising in mathematical physics. The solution procedure of this method, with the aid of Maple, is of utter 

simplicity and this method can easily extended to other kinds of nonlinear evolution equations. 
Ch The present Letter is motivated by the desire to extend the exp-function method to the general types of Burgers-Fisher equation 

with variable coefficients, which reads: 



Ut - Uxx + Oi{t)uUx = I3{t)u{l - u), (1) 

where a{t),j3{t) are arbitrary functions of t. 

The Burgers-Fisher equation has a wide rang of applications in plasma physics, fluid physics, capillary-gravity waves, non- 
linear optics and chemical physics. When a{t) — 0, 13{t) is a arbitrary constant, Eq.([T]l turns to Fisher equation. 



' ut-u.^x = l3u{l - u), (2) 

o ; 

. Kobnogorov, Petrovskii, and Piskunov studied this equation in ||28ll Thev showed if initial datum satisfies some conditions then 
the solution of Eq.® approaches a travelling wave of speed Co = 2. Exact solution of Eq.(|2]) was found by Ablowitz and 
J> Zeppetella in li29il at Co — i"^- When — 0, a{t) is a arbitrary constant, Eq.([T]l turns to Burgers equation. 



ut - Uxx + auux = 0, (3) 

which is used to describe the spread of sound wave in the medium with viscidity and heat exchange if we do not consider the 
medium's frequently dispersive character and the slack comfort process, at the same time, the Burgers equations with variable 
coefficient can be used to describe the cylinder and spherical wave in these questions such as overfall, traffic flow model and so 
on. 

Therefore, there are important theoretic and factual value for us to look for the exact solutions of Eq.([T]l- 



II. THE EXP-FUNCTION METHOD 



We now present briefly the main steps of the Exp-function method that will be applied. A traveling wave transformation 

u — u{£), ^ = kx + wt converts a partial differential equation 

'i'{u,Ut,Ux,Uxx,Utt,- ■ ■) ^0, (4) 
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into an ordinary differential equation 

wu', ku\ k^u", w^u'\ kwu", • • • ) = 0, (5) 
The Exp-function method is based on the assumption that traveling wave solutions can be expressed in the following from 

u{0 = '^""'^''P^"^^ (6) 

where c, d,p and q are positive integers which are unknown to be further determined, a„ and 6,„ are unknown constants. To 
determine the values of c and p, and the value of d and q, we balance the linear term of highest order in Eq. (|5]l with the lowest 
order nonlinear term, respectively. 



III. APPLICATION TO THE BURGERS-FISHER EQUATION 

In order to obtain the solution of Eq. we consider the transformation 

u = u(^), ^ = kx+ / T{t)dt, (7) 



where fc is a constant, T{t)is an integrable function of t to be determined later, then Eq. O becomes an ordinary differential 
equation, which reads 

T{t)u' + ka{t)uu - k'^u" - I3{t)u{l - w) = 0, (8) 

where prime denotes the differential with respect to ^. 

According to the Exp-function method, we assume that the solution of Eq. dS) can be expressed in the form 

^^^^ ^ En=^cQ" exp«) ^ Q_, exp(-ce) + ■ ■ ■ + fld cxp(dg) 
Em=-p ^™ exp(m^) b-p exp(-p^) + h 6g exp(q^) ' 

where c, d, p and q are positive integers which are unknown to be further determined, a„ and are unknown constants. 

In order to determine values of c and p, we balance the linear term of highest order in Eq. (|8]lwith the highest order nonlinear 
term, and the linear term of lowest order in Eq. ^ with the lowest order nonlinear term, respectively. By simple calculation, we 
have 

^ . /^iexp[(d + 3,)g] + ... ^ ^^^^ 
h2 exp(4gC) H 

and 

u.s _ ^3 exp(2d + g)^ + ■ ■ ■ _ fea exp(2rf + 2q)^ + ■■■ 

where hi are determined coefficients only for simplicity. Balancing highest order of Exp-function in Eq. (fTOband ( fTTT l. we have 

d + 3q = 2d + 2q, d = q. (12) 
Similarly to determine values of c and p, we balance the linear term of lowest order in Eq. ([8]) 

."(0 = ■■■ + ;---p[-(;+y , (13) 

■ ■ ■ + S2 exp(— 4p4j 

and 

g3exp[-(2c + p)]e + --- ■■■ + 53 exp[- (2c + 2p)]g 

(0 = ^ r^rT\ — = \ / . c\ — ' (^^^ 

■ ■ ■ + S4 cxp[-6pt,) ■ ■ ■ + S4 exp(-4p4 j 

where Si are determined coefficients only for simplicity. Balancing highest order of Exp-function in Eq. (fTsTland (fT4l) . we have 



c + 3p = 2c + 2p, c = p. 



(15) 
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we can freely choose the values of c and d, but the final solution does not strongly depends upon the choice of values of c and 
d{26\. For simplicity, we set 61 = 1, p = c = 1 and d = q = 1 Eq. ©becomes 

, . _ ai exp(^) + gp + o-i exp(-g) 

exp(0 + 6o + 6-iexp(-0 ' ^ ^ 

Substituting Eq. (fT6]linto Eq. (|8]l, and with the help of Maple, we have 

^[(73 exp(30 + C2 exp(2e) + Ci exp(e) + Co + C_i exp(-e) + C_2 exp(-2e) + C_3 exp(-3e)] = 0, (17) 

where 

A = (cxp(e) + &o + ^'-iCxp(-0)^; 
C3 = -aiPit) ^ aim- 

C2 — —2 Oi6o/3(t) — A:aiaoQ;(t) + a^; bQ/3{t) + ai6oT(i) + /co^6oa(t) — aoT{t) — ao/3{t) + k^aibo — k^a^ + 2 aiao/?(i); 
Ci — aQ/?(i) — 2 kaia-ia{t) — aoboT{t) + 2 aoaibol3{t) — k^aib^ — kaQa{t) + 2kalb-ia{t) — a_i/3(<) — 2a_iT(i); 
+aiblT{t) - 2 aoboP{t) + 2 aia_i;3(i) + fc^ao^o - aibl(i{t) - 4 A;^a_i + 2ai&_ir(t) + 4A:^ai6_i 
-2ai&_i/3(t) + a\b-i(3{t) + kaoaibQa{t); 
Co — 6fc^ao6-i + 2aoa_i/3(t) — 2aibob^i(3{t) — 2a_i6o/3(i) — 3a_i6o''"(i) — 3fc^a_i6o + 3ai&o^-i''"(i) — 3fc^ai6o&-i; 
+2aia2; bo(3{t) + 2aiao&-i/3(t) — 3fcaoa_iQ!(t) + 3fcaiao6_iQ!(i) — 2aofe_i/3(i) — ao6o/3(i) + ao^o/3(0! 
C_i = -2fca^ia(t) - ai6^i/3(t) + a2;^/3(t) + ao6_i/3(t) + 2fca_iai6_ia(t) + 2ai6^iT(t) - a_i6or(t) + 2aQa^ibol3{t) 
-a^ibll3{t) + k'^aob-ibo + 2a_iai6_i/3(i) - 2aob-ibal3{t) - ■ik'^aib\ + aob^ibaT{t) - 2a_i6_i/3(i) 
—kaQa^ibQa{t) — 2a_ife_iT(<) — k^a^ib^ + 4fc^a_i6_i + fcaQ6_ia(t); 
C-2 = -2a-ibQb-il3{t) ~ k^aob\ - aob\(3{t) - ka\boa{t)] +a\bol3{t) + k^a-ibob-i + ao6?.iT(t) 

+fca2ao&-ia(t) — a_i6o&-i'''(i) + 2a_iao62/3(0; 
C_3 = al^b^im-a-ib'-im- 



Equating to zero the coefficients of all powers of exp(n^) yields a set of algebraic equations for ai, oq, a_i, 
&o, bi,k, T{t), a{t), (i{t) . Solving the system of algebraic equations with the aid of Maple, we obtain; 
easel. 

T{t) = -k'^ - P{t), a{t) = ^^, oi = 0, ao = ao, a_i = ao6o, &o = 6-1 = 0. (18) 



ease2. 



ease3. 



T{t) =2k^ + ^l3{t), a{t) = -4k, fli = 1, ao = 0, a_i = 0, bg = 0, 6-1=61. (19) 



T{t) = -k^ - I3{t), a{t)^2k, ai=0, = ^" ^ ^ , a_i = 6_i, 60 = 60, 6_i=6_i. (20) 

ease4. 



T(t)=e+f3{t), a(i) = -2fc, ai = l, = \ a_i = 0, 60 = 60, 6_i = 6_i. (21) 

easeS. 

T{t) = -k'^, a(t)^^Y^, ai = 1, oo = Oq, a_i = -65 + ao6o, 60 = 60, 6_i = 0. (22) 



ease6. 



fc^(12y26l7a;S + 7V26ii + A0V2bt,al + 45a[j6_i + 2al, + 37ao) 
T(i) = 3 5 ; , a[t) = -2k, 6_i = 6_i, (23) 

10V2a26^i + y261i + 6y263;"ao + + 7ao62_i + 15a36_i 
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^^^^ ^ 6fc'(2aobli72 + v/26l7ag + 3ao) 

Substituting Eqs. ([T8]l-(|24]l into Eq. ([T6]l yields 



fli = 1, flo = ao, a_i = 0, 6o = 



(24) 



wi(a:,t) 



exp [/cx - / A;2 + p{t)dt] + bo 



U2ix,t) 



U4(x, t) 



exp 


kx 


+ /2fc2 + ^ 


^dt 




exp 


kx + j2k^ + ^dt 


+ 6_i exp 


—kx 


_ / 2fc2 + mrfi 



(25) 



(26) 



exp [kx - J k'^ + I3{t)dt] +bo + b-i exp [-kx + / fc^ + 



exp [kx + + mdt] + 



U5{x,t) 



ue{x,t) 



exp [fca: + / fc2 + l3{t)dt] +bo + 5_i exp [-/ex - J k'^ + l3{t)dt] 
exp{kx — k'^t) + flo + (— 6q + ao6o) exp(— fcx + fc^t) 



exp(fca; — fc^^) -|- fog 
exp(A;x + / T{t)dt) + oq 



exp(fca; + / T{t)dt) - ^f2hT\ + 6-i exp(-A:x - / r(t)di) ' 



(27) 



(28) 



(29) 



(30) 



where T(t) 



fe^(12^26_iaQ+7\/26^j+40\/26^j^aQ+45aob_i+2ao+37ao 



solutions obtained in LSftl using the first integral method. 



As we will show later, these solution have included all the 



IV. SOME DISCUSSIONS ABOUT THE SOLUTIONS 

(1) If we take &o = in Eq. dZSl l. we have 



u\\(x^ t) = ao exp 



-kx+ / k"^ + (3{t)dt 



(31) 



(2) If we take a{t) = a is a constant, 6_i = 1 and (3{t) — ^"'^ ° in Eq. ( |26] ). where c is a arbitrary constant. Then we have. 



U2i[x,t) - - -tanh 



-(x-rf) 

This is the solution (36) obtained in H30(l . Other solutions in H30(l also included in our solution (25) - (30). 
(3) If we take &o = 4, = 1 and fc = 1 in Eq.dZTli. we have, 

-l + (2 + 2V^) CSC [x — / 1 + Pit)dt] + coth [x- Jl + I3it)dt] 



(32) 



U3i{x,t) 



4csc [x - / 1 + (5{t)dt\ + 2 coth [x - J I + l3{t)dt] 



(33) 



FIG. 1: (a) Solution of Eg.llSTt with ao = 1, fc = 1 and — t. (b) Solution of Ea.l l32l ) with a = 4 and c = 5. 
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FIG. 2: (a) Solution of Eq.lHD with /3(t) = cost. (b) Solution of Eq.l[34ll with I3{t) = -1 + 3sini. 



(4) If we take bo = 0, 6_i = -5 and fc = 1 in Eq. (|28]l. We get, 



U4l{x,t) 



cosh [x + / 1 + + sinh [a; + / 1 + l3{t)dt] - \/5 

-4 cosh [x + Jl + /3{t)dt] + 6 smh [x + J I + /3{t)dt] 



(5) If we take bg = 2 and cq = |, we have 

-"51 = 

(6) If we take oq = in Eq.(l30ll. we have 

U6l{x,t) 



2tanh(fca; - Pt) + | sec(fca; - Pt) 
1 + tanh(A:a:: - Ft) + 2 sec(fcx - k^t) ' 

exp(fca; + 7kH) 



exp(fca; + Tk'^t) — ^26_i + 6_i cxp(— /ct — 7fc^i) ' 
These solutions for local area structure are shown in Figs. 1-3. 



(34) 



(35) 



(36) 



V. CONCLUSION 

The Burgers-Fisher equation with variable coefficients is investigated by Exp-function method. The generalized Travelling 
wave solutions of this equation are obtained with the help of symbolic computation. From these results, we can see that the 
Exp-function method is one of the most effective methods to obtain exact solutions. 

Finally, it is worthwhile to mention that the Exp-function method can be also extended to other nonlinear evolution equations 
with variable coefficients, such as the mKdV equation, the (3H-l)-dimensional Burgers equation, the generalized Zakharov- 
Kuznetsov equation and so on. The Exp-function method is a promising and powerful new method for nonlinear evolution 
equations. This is our task in future work. 
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FIG. 3: (a) Solution of Eq.lElll with fc = 2. (b) Solution of Eq.^ with A: = 1 and 6_i = 2. 
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